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B.P. Leonard* 

Center for Computational Mechanics 
The University of Akron 

SUMMARY 

Judging by errors in the computational-fluid-dynamics literature in recent years, it is not 
generally well understood that (above first-order) there are significant differences in spatial 
truncation error between formulations of convection involving a finite-difference approxima- 
tion of the first derivative, on the one hand, and a finite-volume model of flux differences 
across a control-volume cell, on the other. The difference between the two formulations 
involves a second-order truncation-error term (proportional to the third-derivative of the 
convected variable). Hence, for example, a third (or higher) order finite-difference approxi- 
mation for the first-derivative convection term is only second-order accurate when written 
in conservative control-volume form as a finite-volume formulation, and vice versa. 

FINITE-DIFFERENCE AND FINITE- VOLUME FORMULATIONS 

Consider the model constant-coefficient one-dimensional pure convection equation for a 
scalar <t> 

*4 * u - S(x, ,) (1) 

dt dx 

where 5 is a known source term, and assume that a numerical solution is sought using a 
discrete grid of constant step-width h. As usual, let </>, represent the numerical approxi- 
mation of <£ at grid-point i. 

*Work funded under Space Act Agreement NCC3-233. 
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( 2 ) 


A finite-difference formulation of Equation (1) attempts to simulate 



and, in particular, the spatial first-derivative convection term is written in terms of node- 
values of <£. The modelled first derivative is then equal to the true first derivative at i, plus 
truncation error terms: 


H\ 

dxl 


model 



( 3 ) 


The leading term in (T.E.) FD (i.e., the term involving the lowest power of h ) is conven- 
tionally called the “order” of the finite-difference discretization. 

On the other hand, consider integrating Equation (2) with respect to x, from -h/2 to 
+ h/2, and dividing by h. This gives 


H t 

~dT 


u (<l> r - <f> t ) 

h 


+ Sff) 


( 4 ) 


where the bars refer to spatial averages, and left and right control-volume face-values are 
indicated. This is the finite-volume formulation of Equation (1). 

In this case, one writes 


H* r $ 1 ) model 

h 


Hr ~ <M 

h 


+ (T.E.) pv 


( 5 ) 


where the right-hand side involves the true face- value difference. Once again, the leading 
term in (T.E.) PV is the order of the finite-volume discretization. 
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It is often assumed (especially in recent CFD literature) that, if a finite-difference model is 
written in flux-difference form, then (T.E.)fd is the same as (T.E.) PV . But, as will be 
shown, except for the leading term in first-order formulations, 

(T.E.) fd * (T.E.) FV (6) 

The confusion is apparently based on the fact that the finite-difference model of the first 
derivative can often be split into two parts; i.e., 

I d c/> \ = jy - / (7) 

\ 3x / model h 

where 1), and this is sometimes treated as a finite-volume formulation (with 

the assumption that the truncation error is the same). But if Equation (7) is to be treated as 
a finite-volume model, one must recompute the truncation error according to Equation (5). 

FACE-CENTERED TAYLOR EXPANSIONS 

For definiteness, consider the classical second-order central finite-difference approximation 
for the first derivative: 


/ d<M = Ki ~ (8) 

\ dxl model 2 h 


First, make Taylor expansions about grid-point i. For example, 



</>,♦ i = <t> t + <t>i h + ± <t> 

0t | Ml * 

i hr + - fa h s + ... 

(9) 

and 





4>i - 1 = <t>i ~ 4>i h + i <t> 

t h 7 - 1 4.7 h> + ... 

(10) 

so that 





<A 1+ i “ ^«-i = 2< t>i h + } <; 


(11) 
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thus giving the well-known result that 


fy+i ~ 4 > t - 1 

2h 


= (4*) ♦ 1 * -l *, w h' 

\dxli 6 120 ' 


( 12 ) 


verifying that this is, indeed, a second-order approximation to the first derivative. 
But this model can be rewritten in the form of Equation (7) by identifying 


and 





<f>i + <t>,. i 
2 


(13) 


(14) 


In other words, the modelled left and right face-values are taken to be just the arithmetic 
means of the node-values on adjacent sides of the individual faces. Note, as required by 
conservation, that 


= *;(/- 1) (15) 

Now the model can be considered as a finite-volume formulation simply by writing 


l^t * i + + 4>t- 1\ 

*'r _ V 2 I V 2 I (<t> r ~4> t ) ^ c , 

h h h <T - E - ) "'' 


(16) 


In order to assess the truncation error, expand the node-values about individual control- 
volume face locations: 


<t> 


i* 1 



(17) 
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( 20 ) 


( 21 ) 


4>i * <£ M 


i 4>; h 2 + _L ^ iv> /i 4 

8 1 384 1 


so that 

<.*: - *;> _ <*, - *,) ^ i /«; - *;'i ^ i i4 (23) 

ft A 8 l ft / 384 \ ft / 

But, from Equations (17) and (18), 


• 4 2 h y*‘ ) ♦ 


and similarly for 4>'i . Then, using Equation (9) together with the following expansions of 
face- values about grid-point i, 


5 



( 25 ) 


and 



the difference of face-second-derivatives appearing in Equation (23) can be written as 



h 


4>i 


i 

24 




— 
1920 


(vii) 


(27) 


giving 


+ 4>i\ l<t>i + </>,-! 


(*, - <M , 

- 1 * T *> 


1 

128 


*!" h‘ 


(28) 


Thus, by comparing Equation (12) and (28), one sees that 


(T.E.) pd = i *«** * ... (29) 

whereas 

(t.e.) fv = * ^ . (30) 

This, of course, is a significant difference, even though both formulations are second-order 
accurate. Note that the difference in the truncation errors is 

(T.E.) fd - (T.E.) pv = JL t'i'h* * JL h* * ... (31) 

and a result similar to this will be found in general to be true for any convection formula that 
can be simultaneously viewed either as a finite-difference formula for (30/dx), or a finite- 
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volume formula for (4> r - <f> t )/h. In fact, referring to Equations (25) and (26), continued 
through fifth-order, one finds that, irrespective of the numerical scheme, 



which explains the difference between Equations (29) and (30). 


OTHER COMMON DISCRETIZATIONS 

In addition to the second-order central-difference formulation considered above, it is 
convenient to summarize a number of other discretizations commonly used in convective 
modelling. 

First-Order Upwinding 

For u > 0, the convective term in Equation (2) is written 


- u 




(33) 


From Equation (10), viewed as a finite-difference formulation, this gives 


(ti ~ <K_ i) 

h 




(34) 


which, as expected, is first-order accurate. Viewed as a finite-volume model, for u > 0, 
the face-values are written (with upwind bias) as 


(^r^model 



(35) 
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and 


(^f)n nodel ^/-l = $1 (^) + ••• 

And this gives 

(f “ 'Mmodd (<^r “ 4,) 1 (f ~ f\ , 

A = A " 2 \ A ) h + 


(36) 


(37) 


But, from Equations (17)-(20), 


4>r = 


_ (^1*1 - <*>,) 1 J.”' fc2 

r - 17 *r h 


and 


so that 


& _ (*^< ^i-l) 1 A '" h 2 A 

4> t - j + 


(38) 


(39) 


<t>'r - 4>'t _ (* M - 2<t> t + 0 M ) , (f - f) 


24 


h 2 + 


(40) 


and the second central-difference can be written 


(*,♦. " 2 f + 0 M ) 

h 2 


f, * ± A 1 + 


as is well known. This means that Equation (37) becomes 

(4»r - )modd (f ~ <l>,) l u 

r - T " X Pi h + 


(41) 


(42) 


so that the leading truncation error is the same as that of the finite-difference formula, 
Equation (34). This, of course, is to be expected from Equation (32). 
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Second-Order Upwinding 

For u > 0, if one interpolates a fully upwind-biassed parabola through i, (/ — 1 ) , and (i-2), 
the corresponding first-derivative at / is 

/a^\ = (H - jjV-i + 4> t - 2 ) (43) 

\ (3X / model 2 H 

. (W) + (44) 


But the right-hand side of Equation (43) can also be written in finite-volume form as 



(* r - <f>t) 

h 



\ ^ 
4 


(iv) 


(45) 


which again conforms with Equation (32). Note that, in this case, face- values are obtained 
by linear extrapolation from upwind nodes. 

Third-Order Upwinding 

This time, for u > 0, interpolate a (partially upwinded) cubic through (i + 1), i, (i-1), and 
(i -2) . The corresponding first-derivative at i is then 

/ d<M = (<I> M ~ _ (Ki ~ ty- H-i ~ jVg) (46) 

\ dxl model 2 h 6 h 

Written in this form, one can see that the third-difference will cancel the leading truncation 
error in Equation (12), giving 

(£L ■ *:•’*<*... <«> 

which is indeed a third-order accurate representation of the first-derivative at i. 
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On the other hand, Equation (46) can be rewritten in finite- volume form by identifying face- 
values (for u > 0) as 


(Vrfm odd ~ ~ ~ T W\ + l " 2< £< + ^i-l) 


(48) 


and 


(J. s _ (4, + 0 <-,) 1 f. ^ . 

W*/) model ” S “ T (<£, ~ + 0,.. 2 ) 


(49) 


But this gives 


($ r ^f)model 


(^V ) 1 j. j,2 1 ^(»v) ^ 3 J 

h ~ Ta ^ h + 71 +* k + 


(50) 


which of course is only a second-order accurate approximation. 

In order to achieve a third-order accurate finite-volume representation, one needs to annihi- 
late the leading truncation error in Equation (28). This is achieved by writing (for u > 0) 


(*,). 


model 


— — " I tot* i “ 2 ^« + 4>i- 1) 


(51) 


and 


(0 


f) model = — — tT~^~ ~ 7 (<t>i ~ 2 ^f-l + ^i- 2 ) 


<52) 


giving 


(*r - M 


model 




-L 4 
16 ‘ 


(hr) 


— 

12 * ' 


(v) 


(53) 


which is seen to be third-order accurate. Equations (51) and (52) represent the well-known 
QUICK formulas for face-values, obtained by interpolating a parabola through the two 
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nearest node-values together with that of the next adjacent upwind node. In summary, the 
1/8 factor on the second-difference terms is appropriate for a finite-volume formulation, 
whereas the 1/6 factor corresponds to the finite-difference model of the derivative. In 
practice, the difference between using 1/8 and 1/6 (in a finite-volume formulation) is 

observed to be quite small. Note that second-order upwinding can also be written in a 

similar form, using a factor of 1/2 on the second-difference terms. In this case, however, 
results are significantly less accurate. 

Higher-Order Formulations 

The simplest way to construct higher-order formulas is to start with a known formula and 
add higher-order difference terms to cancel the leading truncation error. For example, if one 
were trying to construct a fourth-order accurate approximation to the first-derivative at i, the 
appropriate formula would cancel the </>, term in Equation (12) without introducing an /i 3 
term. This can be done by using the average third-difference centered at node i given by 

[(+1.2 - 3+1.1 + 3+1 - +,_,) * - 3 < t >, * 3+i_, - +,_ 2 )] 

- -i (+i.2 - 2+,., * 2+1-, - +,-,) (54) 


so that 


l d<f> \ m (4 > m ~ ^j-i) _ (</>„2 ~ 1 + ~ jg) 

\ dxl model 2 h 12 h 


(55) 


On the other hand, the appropriate fourth-order finite- volume formulation would use the 
average second-difference centered at a face. For example, 


\_ 

2 


[(<£,•♦ 2 " 1 


</>/) + (4>i . i ~ 2( t>i + 4> /_i)] 


^ Wi* 2 ~ + ^i-l) ( 56 ) 
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so that the appropriate face- value is 


(^r) model 


1 + ^i) _ (^< + 2 “ $1+1 ~ 4>i + ^i-l) 
2 16 


(57) 


with a similar formula for <j>, (reducing all indexes by 1). 

Once again, one sees that Equation (55) could be rewritten in finite- volume form using 


(<M 


model 


(^<♦1 + 4>i) _ i+2 " i+1 ~ </>, + ^,-i) 

2 12 


(58) 


with a similar formula for the left face. But this would result in a finite-volume formulation 
that is only second - order accurate, as predicted by Equation (32). 


CONCLUSION 


Equation (32) shows that there is a significant difference between the first derivative at a 
node and the face-value difference (divided by h) across a control-volume cell. If a 
convection scheme is constructed on the basis of modelling ( d<f>Idx) i , with truncation error 
(T.E.) fd , and then rewritten in conservative finite-volume form, the truncation error must 
be recomputed according to Equation (5), using Taylor expansions about face values. The 
difference in accuracy shows up in steady-state calculations, where dfy/dt = d4>Jdt = 0. 
Interestingly enough, if one writes, in the vicinity of grid-point i, 


<K*) = <f>. + 4>'i x 


i * f 

+ \4> t x 2 


1 A," 
6 +* 


-L 0f V) + 


24 


(59) 


and then computes the control-volume cell average 



the result is 


( 61 ) 


This, for example, explains the difference between the 1/8 factor in the third-order steady- 
state QUICK scheme and the 1/6 factor in the third-order time-accurate QUICKEST scheme, 
which was pointed out thirteen years ago 1 . 
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